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Abstract A review of the low temperature properties of Kondo lattice systems is 
presented within the mean-field approximation, focusing on the different charac- 
teristic energy scales. The Kondo temperature, Tk, and the Fermi liquid coherence 
energy, 7b, are analyzed as functions of the electronic filling, the shape of the non- 
interacting density of states, and the concentration of magnetic moments. These 
two scales can vanish, corresponding to a breakdown of the Kondo effect when an 
external magnetic field is applied. The Kondo breakdown can also be reached by 
adding a superexchange term to the Kondo lattice model, which mimics the intersite 
magnetic correlations neglected at the mean-field level. 



1 Introduction 

Rare-earth and actinide based compounds exhibit extremely rich phase diagrams, 
with signatures of heavy fermion behavior, unconventional magnetism, or super- 
conductivity (HE). At high temperature, the main physical properties of these sys- 
tems are well reproduced by single impurity models, which describe the coupling 
between conduction electrons and one 4/ or 5/ ion. For dense systems, the single 
impurity models fail to describe the low temperature properties, which are charac- 
terized by the formation of a non local coherent state. In this regime, models with a 
periodic lattice of / ions are more appropriate. 

Here, we consider more specifically dense compounds where the / orbital is oc- 
cupied by one electron (Cerium) or one hole (Ytterbium). In the low temperature 
regime where the crystal field splitting lifts the degeneracy of the / orbital, these 
impurities are modeled by effective local spins Si = 1/2. The system is thus de- 
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scribed by the Kondo lattice Hamiltonian, 

tf = E( £ k-M)cL c ka+^I> S '' (1) 
kc ! 

where c^ a (ck a ) describe creation (annihilation) operators of conduction electrons 
with spin a =|, J, and momentum k. The Kondo interaction results from a local 
antiferromagnetic coupling Jk between the density of spin of conduction electrons at 
site i, Sj, and the Kondo impurities, S;. The chemical potential jU fixes the electronic 
filling to n c conduction electrons per site. 
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Fig. 1 Schematic phase diagram. The reduced temperature T/Tk is plotted versus Tq/Tk, which 
is considered a tunable parameter that can vary with the electronic filling, impurity concentration, 
magnetic field, or the shape of the non-interacting DOS. Figure from Ref. |4] 



The Kondo model has been extensively studied throughout the last decades |3). 
At high temperature, the Kondo interaction can be considered a small perturbation: 
conduction electrons and Kondo ions are weakly coupled. The transport proper- 
ties, which are determined by the conduction band, correspond to those of a normal 
metal. The magnetic susceptibility, governed by the Kondo free moments, has a 
Curie-Weiss form. The entropy is large, of the order of In 2 per site. A crossover 
occurs at the Kondo temperature, Tk, below which the Kondo interaction cannot 
be treated by perturbative methods. Experimental signatures of this crossover in- 
clude, for example, a logarithmic increase of the resistivity when the temperature 
decreases. Other signatures involve a saturation of the magnetic susceptibility, and 
a significant decrease of the entropy. At lower temperature, if we neglect magnetic 
or superconducting instabilities, the physical properties are characteristic of a uni- 
versal heavy Fermi liquid: the specific heat vanishes linearly with the temperature, 
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Cy{T) w yT; the local magnetic susceptibility, %{T), as well as the resistivity are 
constant at T = 0, with quadratic (i.e., T 2 like) variations at low T. This low tem- 
perature regime is characterized by an energy scale, the coherence temperature, 7b, 
which can be equivalently determined from the specific heat Sommerfeld coeffi- 
cient, To = l/y, or the zero temperature magnetic susceptibility, Tq = l/x(T = 0). 
Figure Udepicts the schematic phase diagram of the model, as a function of Tq/Tk- 
The determination of 7b and Tk from the entropy is illustrated by Fig.|2] The connec- 
tion between 7b and the thermal and electric transport properties have been analyzed 
byZlatic etal. 016). 

Note that the authors of Ref. Q discussed two energy scales: one of them corre- 
sponds to our definition of the Kondo lattice temperature, Tk, from the temperature 
dependence of the entropy [see Fig. |2). The second is the single impurity Kondo 
temperature, which, within the mean-field approach, is equal to the lattice Kondo 
temperature. The coherence temperature 7b that we analyze here is not considered 
in Ref. Q. 



Se-(T) 




Fig. 2 (Color online) Schematic plot of the electronic contribution to the entropy S e - as a func- 
tion of the normalized temperature T /Tk for three cases: Tq» Tk (red dash dotted line), Tq rs Tk 
(blue solid line), and Tq «Tk (green dashed line). The dotted lines indicate the linear Fermi liquid 
regime S e - (T) = T /Tq, with a rescaled slope Tk/Tq. For T > Tk the three curves are identical, re- 
flecting the linear contribution from the conduction band S e - (T) = ln2+ T/D, The black dot refers 
to a standard experimental determination of Tk from the electronic entropy: S e - (Tk) = pln2. On 
this schematic plot, we used p=l, which coincides with the vanishing of the effective hybridiza- 
tion, r = 0, denning Tk within the mean-field approach. Experimentally, where Tk is a crossover, 
one can choose, e.g., p = 1/2. Figure from Ref. |4]. 



The physical properties of systems with a small concentration of magnetic ions 
are universal and characterized by a single energy scale, 7b = Tk. The identity re- 
lating the coherence and the Kondo temperatures is consistent with the exact solu- 
tion of the single impurity model. The situation is different for dense systems for 
which more than one energy scale can be identified. As an example, the Ruderman- 
Kittel-Kasuya-Yosida (RKKY) interaction, Jrkky, can compete with Tk, and a mag- 
netic instability can be obtained for a lattice system j8j. For a long time, it has 
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been believed that 7b and Tk remain equal to each other in the heavy Fermi liquid 
phase of a Kondo lattice model. In this context and in response to discrepancies be- 
tween theory and photoemission spectroscopy on rare-earth-based compounds with 
respect to the scaling with Tk, alternative scenarios involving phonons have been 
proposed (9j. On this basis, whether the Kondo model is capable of predicting pho- 
toemission spectroscopy on YbAl3 has been a subject of controversy iflOl fTTl PPUl . 
The specificity of Kondo lattice systems, with 7b different from Tk, has been con- 
firmed by experiments including magnetic susceptibility, specific heat, Hall coeffi- 
cient measurements and X-ray absorption spectroscopy in Cerium, Cei-.vLajJ^Ge?, 
CeIr 2 - A (Rh,Pt) I Ge 2 , CeIr 2 Ge 2 -^(Si,Sn). v HD, CeNiSi 2 O, and Ytterbium com- 
pounds, YbXCu 4 (X=Ag, Cd, In, Mg, Tl, Zn) 02), Ybi-.^Lu^Als lfl6l . 

The first suggestion that 7b could be much smaller than Tk was discussed by 
Nozieres in the framework of the exhaustion problem irT71[T8l . This prediction stim- 
ulated complementary theoretical works, including approximated methods based 
on a strong coupling approach |fl9l or mean-field calculations l20l |2T1 l22l [23), 
and numerical simulations 11241 l25l l26l l27l |2§1 using dynamical mean-field theory 
(DMFT) [29 30). Even if the initial prediction of Nozieres turned to be quantita- 
tively wrong IF2TII3TI . all the theoretical calculations converge to the same qualitative 
conclusion: in a Kondo lattice, 7b can be different from Tk, and the ratio between 
these two energy scales can be tuned, for example, by varying the electronic filling 
of the system. Finally, these theoretical works shed some new light on the formerly 
controversial photoemission analysis. 

The aim of this work is, first, to review how the two energy scales, 7b and Tk 
depend on physical parameters: electronic filling, shape of the density of states, 
concentration of magnetic ions. This is done within the simplest relevant approxi- 
mation for the Kondo interaction: a mean-field decoupling. Then, we review how, 
within the mean-field approach, the Kondo phase can be destabilized by a magnetic 
field or by magnetic inter-ion interactions. 



2 Mean-field formalism 

The mean-field approximation for the Kondo lattice was first introduced by Lacroix 
and Cyrot (32). It was reformulated by Coleman lf33l and by Read, Newns and 
Doniach 11341 . as a large— N approximation for the N— fold degenerate Coqblin- 
Schrieffer model ll35l . It was shown that magnetic instabilities require expansions 
up to the order 1 /N, i.e. fluctuations around the mean-field. Nevertheless, the heavy 
Fermi liquid phase is well described in the limit N = °°. The analysis of To and Tk 
can thus be already performed at the mean-field level. Here, we describe the main 
lines of the mean-field approximation for the Kondo lattice Hamiltonian ([T). 

The Kondo impurities are represented by local auxiliary fermions as follows: 
Sf = \{f^fi] - fl\fil)> s l~ = fj\fib and s i = fJifi\- The Kondo interaction is thus 
rewritten as JkSiSi i— > Y,a& c 'io c '<T'flo>f icr > which describes the spin-flip processes 
between conduction electrons and local moments. This mapping is exact as far 
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as the Hilbert space is restricted to the sector of one auxiliary fermion per site, 
fr\fi1 +fJifi\ = !• The mean-field solution is obtained within the two following ap- 
proximations: (i) The local occupation of the auxiliary fermions is equal to one only 
on average. This corresponds to describing the Kondo spins by an effective local 
/—level that is half full. This effective filling is controlled by introducing a second 
chemical potential, A. (ii) The Kondo interaction is replaced by an effective one- 
body term, obtained from a mean-field decoupling of the f and j. components. The 
mean-field approximation for the Kondo lattice Hamiltonian ([TJ yields 

H = E( £ k -^claCkc + rYjclfia+f^Cia] " , (2) 

kc ia ia 

where the effective hybridization is determined by the self-consistent relation 

r =£-p (3) 

Here jV is the number of lattice sites. The chemical potentials for c— electrons, 
ju, and /—fermions, A, are determined by the constraints n c = -\Y*ia{ c ]o c io) an d 
1 = -ipYtioifiafict)' The mean-field Hamiltonian (ffj) describes an effective system 
of conduction electrons hybridized with local / levels. The correlation effects are 
renormalized into the self-consistent hybridization r. This approximation captures 
two important features of the Kondo lattice: at high temperature, we find r = 0, and 
the system is described as a paramagnetic light metal (c electrons) decoupled from 
local free moments (/—fermions). This picture is oversimplified but it succeeds in 
describing qualitatively the experimental situation, in which, above the Kondo tem- 
perature, conduction electrons are weakly coupled to the local moments. Within the 
mean-field approach, the Kondo temperature Tg is thus defined as the temperature 
for which a hybridization r / occurs. At very low temperature, the physical proper- 
ties correspond to a Fermi liquid and the excitations of the system correspond to the 
creation of non-interacting, heavy, fermionic quasiparticles. The latter are a linear 
combination of the light c— electrons and the heavy /—fermions. The Fermi-liquid 
regime is characterized by an energy scale, 7b, which can be defined identically 
from different physical properties of the ground state: the quasiparticle density of 
states, p = 1 /7b, the Sommerfeld coefficient, F = limr^,oCv{T)/T = 1 /7b, or the 
local impurity spin susceptibility, Xloc(T =0) = 1 /7b. We have obtained explicit 
expressions for 7g and To in the limit of small Jg 141 I2TI : 

T K = F K [n c ,po\e- l l jKpQ{iMs) , (4) 

and 

T =Fo[n c ,p }e-V J xP<>^) ; (5) 

which depend on the Kondo coupling Jg only within the non-analytic exponential 
factor e _1 / J tPo(w>), The prefactors Fg and Fq are functions of the electronic filling, 
n c , and the non-interacting density of states (DOS), Po((o) = Lk^fw — £k)- Here, 



6 



S. Burdin 



Ho is the non-interacting chemical potential, corresponding to an electronic filling 
n c = ft%Po(G>)d6). For the sake of clarity, we do not write the explicit expressions 
of Fk and Fq here, as they are given in Refs. J4]|2l]. The exponential factor in both 
Tq and Tk, results in a very sensitive dependence of these energy scales with respect 
to small changes in the system (e.g., pressure or doping). Since the ratio Tq/Tk does 
not depend on Jk, it is more likely to be analyzed experimentally as a universal 
function. The mean-field solution yields 2) 



7b 



D + Mq \ 1/2 F shape 
D-Ho) ap (n L )An 



(6) 



Here, D is the half -bandwidth of the non-interacting DOS. The interacting chemical 
potential, corresponds to a filling n c + 1, and the energy Ajj, = — jj$ is related 
to the enlargement of the Fermi surface, a = 1.13 is a number. F s hape is an explicit 
function depending on the shape of po as follows: 



Fshape — exp 



<D+IM 



_i 1°-^ \ po(Mo + <»)-po(Mo) i(3 ' 

: ) 2J-( D+lH)) ) \co\po{ll ) 



(7) 



3 Tuning 7b and Tk 



3.1 Variation of electronic filling 

The effect of electronic filling on the temperature scales of the Kondo lattice has 
been discussed by Nozieres who suggested a possible exhaustion problem |fT71[T8ll . 
Considering N$ Kondo spins coupled to N c < N$ conduction electrons, Nozieres 
started with the following remark: the magnetic entropy which can be released at 
the temperature T « Tk by the formation of incoherent Kondo singlets (i.e., by the 
single impurity Kondo effect) is AS = N c In 2. However, the formation of a coherent 
Fermi liquid ground state is characterized by a vanishing entropy, S(T) s» T /To. The 
freezing of the remaining entropy, S(Tk) ~ (Ns — N c )]n2, thus requires a collective 
mechanism and can lead to To <C Tg when N c <C Ns- 

The exhaustion problem remained for several years an open issue, and its first so- 
lution was obtained within a mean-field approach Ell 1311 . which provides a quan- 
titative description of the electronic filling effect. First, from the analytical expres- 
sions (@]|5j, Tk and 7b have the same Jk dependence, e _1 MtPo(rt)) factor, and the 
ratio Tk/Tq does not depend on the Kondo coupling. This is in contradiction with 
the result of Nozieres, who predicted a ratio T /T K « T K /D « g-V/irPoCtt)) fT7lfT8l . 
The mean-field result, Tq/Tk independent of Jk, was later confirmed by DMFT cal- 
culations [36), and finally accepted by Nozieres 1311 . Nevertheless, Tq and Tk define 
two energy scales with different dependencies with respect to the electronic filling, 
as depicted by Fig. [3] where 7b <C Tk in the limit n c — > 0. The mean-field result 
obtained here is remarkably similar to the one obtained within DMFT combined 
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Fig. 3 (Color online) Tk/D (black solid line) and Tq/D (red dashed line) versus electronic filling, 
for the Kondo lattice. Numerical result obtained within the mean-field approximation, for a semi- 
elliptic non-interacting DOS and Jk/D = 0.75 1211 . 



with Quantum Monte Carlo simulation [see Fig. 1 in Ref. J26]]. The filling effects 
can also be analyzed from the expression ©. Neglecting the band shape effects 



[discussed in section [3721 . we find Tq/Tk p 
approaches the band edge — D, i.e., when n c 



0. 



1/2 



which vanishes when /io 
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Fig. 4 (Color online) Electronic occupation ;v(£k) f° r T/Tg =1.0 (black dashed line), T/Tk = 0.5 
(black dotted line), and T /Tk = 0.005 (red solid line). Numerical result obtained within the mean- 
field approximation for a semi-elliptic non-interacting DOS, n c = 0.15, and Jk/D = 0.75. jUo and 
indicate the chemical potential corresponding to a small and a large Fermi surface, respectively. 
Inset: focus around fl^. Figure from Ref. 1211 . 
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The analysis of the electronic occupation, n c (£k) = (ctck), provides an important 
insight for understanding the physical mechanism leading to two different energy 
scales. For T « Tk, the mean-field result, depicted by Fig. |4] looks like a Fermi 
distribution with a thermal window around the non-interacting chemical potential, 
flo, corresponding to n c electrons (small Fermi surface). For T <C Tk, in the Fermi 
liquid regime, the distribution is spread and forms a step around ji^ (large Fermi sur- 
face). Figure |4]only describes the occupation of the c— electrons, which is fixed to 
n c . At T = 0, there are states of given momentum k which are not fully occupied by 
c— electrons. The Fermi liquid picture is recovered because the quasiparticles are not 
pure c— states, but a linear combination of c and /—states. The quasiparticle occupa- 
tion is complete, i.e., equal to one, for states with an energy £t < fii, and it vanishes 
for states with higher energy. This behavior is consistent with the Luttinger theo- 
rem which predicts that, at T = 0, the Fermi surface contains both c— electrons and 
/— fermions. It is not surprising that the Kondo lattice satisfies the Luttinger theo- 
rem within the mean-field approximation. The reason is that the effective mean-field 
model (f2j) does not contain an explicit many body interaction term. We expect this 
result to survive beyond the mean-field in the Fermi liquid phase. The enlargement 
of the Fermi surface might be a key point in the origin of the difference between 7b 
and Tk- the Kondo temperature is associated with the incoherent scattering of the 
conduction electrons which are in the Fermi window of width Tk around jUrj [see 
Fig- El ■ The coherence temperature, 7b, results from the Kondo effect, but, unlike 
Tk, characterizes a Fermi liquid with a large Fermi surface. 



3.2 Shape of the non-interacting density of states 




Fig. 5 (Color online) Schematic plot of the non-interacting DOS. (a) For a regular DOS and far 
from the electronic half-filling Tq <C Tic. Close to the half-filling, the shape of the DOS around 
/! is crucial for determining Tq/Tk'- (b) po(co) is nearly constant around ft) = fi and To ~ Tg. 
(c) jX is close to a minimum of po and T Q 3> T K . Here, fig « ix indicates the chemical potential 
corresponding to n c non-interacting c— electrons (small Fermi surface), and \i L is the chemical 
potential corresponding to n c + 1 non-interacting c— electrons (large Fermi surface). Figure from 
Ref. |U. 

Here we consider the effects due to the shape variation of the non-interacting 
DOS, Po(®) ID. In order to separate this effect from the electronic filling effects, 
we assume that n c is close to 1, but not exactly half full, so that the system is metallic. 
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Eqs (|6]l and (0 are thus simplified to 




F shape ~ exp 



" f D -^ Pq{ho + (q)-Pq(iM)) 
J-iD+to) 2|o|p (jUo) 



(8) 



A constant po gives Tq ~ Tk, which explains the T /Tk scaling observed in some 
compounds. If )1q is close to a local maximum of po(co), the integrand in Eq. (fSJ 
is negative in the main part of the integration range, and Tq <C Tk- In the opposite 
situation, when /1q is close to a local minimum [see Fig. 0(c)], we find Tq 3> Tk, 
which can be understood by the following argument. The incoherent Kondo screen- 
ing which begins at T w Tk involves a small number of conduction electrons which 
are in the Fermi thermal window of width Tk, around /1q. At lower temperature, 
T < Tk, the Fermi surface is enlarged, due to the contribution of the /—electrons 
[see Fig.Hl. This results from a non-zero hybridization, r ^ 0, in the Kondo phase. 
The formation of a coherent Fermi liquid ground state thus involves all the states of 
the large Fermi surface. In the situation described by Fig.[5](c), where /lo is close 
to a local minimum of the DOS, the further we are from Hq, the more states are 
available for the formation of the coherent Fermi liquid. In this case, there is a kind 
of self-amplification resulting in Tq 3> Tk- With this picture, the magnetic screening 
of the Kondo impurities involves not only the c— conduction electrons but also, in a 
dynamical way, "growing" quasiparticles which contain some /—components. 



3.3 Substitution of magnetic ions 

Here, our analysis focuses on the concentration of Kondo impurities, x. In rare-earth 
based Kondo systems, the latter can be varied, e.g., by Ce-La or Yb-Lu substitution. 
For this purpose, we have introduced a Kondo alloy model, which is a generalization 
of the Hamiltonian ([U where each site of the lattice can randomly either contain 
a Kondo impurity, with a probability x, or not ll37l l38l . The Kondo interaction is 
then treated within the mean-field approximation, and the different configurations 
of impurity distributions are averaged using a generalization of the coherent poten- 
tial approximation [39 1. Results described here remain valid beyond the coherent 
potential approximation [38). 

Figure |6]illustrates the evolution of Tq/Tk as a function of the impurity concen- 
tration, x, for different values of electronic filling, n c . The dilute limit of the model, 
x <C 1 , reproduces the universal behavior of the single impurity model, with Tq = Tk- 
The crossover to the dense Kondo-lattice regime, with Tq ^ Tk occurs for x rj n c . 

The data presented in Fig.|7]are identical to those in Fig. [6] with a rescaled axis, 
x — > n c /x. From this rescaling, an exhaustion regime is identified, for n c < x, where 
Tq/Tk is a universal function of n c Jx. Note that this numerical result has been ob- 
tained here for a non-interacting semi-elliptic DOS which mimics, at low n c , the 
band edge of a three-dimensional system. We expect the universal exhaustion regime 
to depend on the dimension via the exponent characterizing the vanishing of the 
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Fig. 6 (Color online) Tq/Tk as a function of the impurity concentration. From top to bottom, 
n c = 1;0.8;0.6;0.4;0.2. The solid lines are included as a guide for the eye. The curves have been 
computed for a semi-elliptic non-interacting DOS and Jk/D = 0.75. From Ref. 1371 . 



Non universal regime 




Fig. 7 (Color online) Same data as Figure [6] with a rescaled horizontal axis which represents 
n c /x here. Each curve corresponds to a given electronic filling, n c , as indicated by the legend. 
The parameter n c /x has been tuned by varying x, and the curves are cut-off by the finite physical 
minimal value n c /x = n c . 



DOS at the band edge. Apart from the possibility of magnetic ordering, which is 
not considered here, the exhaustion regime might be difficult to access experimen- 
tally for the following reason: the maximal concentration of Kondo impurities is 
x = 1. The only way to reach n c <C x thus involves decreasing the electronic filling. 
Yet, in rare-earth compounds, conduction electrons usually involve more than one 
band. With a single band description, the exhaustion regime would be observable 
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by increasing the concentration of magnetic ions, x. A multi-band system, however, 
allows for another scenario to materialize as x increases: instead of exhausting one 
conduction band, the system might energetically prefer transferring electrons from 
other bands, preventing the filling from accessing the exhaustion regime, n c <C x. 
In this case, a single band Kondo lattice model would not be appropriate anymore. 
Nevertheless, we expect that other Kondo systems can be realized experimentally, 
in which the universal regime n c <C x could be observed. For most of dense com- 
pounds, x « 1, the electronic filling is n c w x. This is a non-universal regime, where 
all the lattice structure becomes relevant. 



4 Destabilizing the Kondo phase 

Microscopically, the Kondo effect is characterized by the formation of local sin- 
glets. The Kondo phase can be destabilized by an external magnetic field, or by the 
fluctuations of the internal Weiss field induced by RKKY interactions. Here, these 
two situations are analyzed from the mean-field approximation: the breakdown of 
the Kondo effect is identified to a continuous vanishing of the / — c effective hy- 
bridization, r, which occurs at T = 0. 



4.1 Effect of a magnetic field 

We consider an external magnetic field, h, applied to the Kondo spins in the longi- 
tudinal direction z. This mimics a situation where the Lande factor of the magnetic 
ions is much larger than the one of c— electrons. A supplementary contribution is 
thus added to the Hamiltonian (Q]i, H ^ H — h^Sf. The mean-field approximation 
described in Sec.|2]is generalized J4), resulting in an effective Hamiltonian formally 
similar to Eq. (0, with spin-dependent /— fermion potentials, A i— > A + ah. The 
effective hybridization, r, remains spin-independent. 

For a sufficiently small magnetic field, a low temperature solution with r ^ is 
obtained. This situation, depicted by Fig. [8] characterizes a phase where the Kondo 
effect coexists with a partial polarization of both Kondo spins and c— electrons. 
A finite critical field h c (T) is obtained, above which c— electrons decouple from 
fully polarized local moments. At T = 0, in the weak coupling limit, Jg <C D, the 
mean-field approximation yields the universal relation = h c (T = 0) = T^/a. 
Since a = 1.13, we have h®. m Tk- This result is not surprising if we consider 
the Kondo temperature as the energy scale characterizing the local singlet for- 
mation: the Kondo effect is destroyed when the Zeeman energy becomes larger 
than Tk- At finite temperature, assuming a constant c— DOS yields the critical line 
[h c (T)/h^ + [T/T K ] 2 = l. 

The T = magnetization m z (h) = ^Ydiifnfi]) ~ {f][fi[)) obtained from the 
mean-field solution is plotted in Fig. [8] as a function of the reduced magnetic field 



12 



S. Burdin 



m z (h) 

Light metal (c electrons) 




Fig. 8 (Color online) Left: Schematic phase diagram of the Kondo lattice as a function of a mag- 
netic field h. The red solid line indicates h c (T) which separates the Kondo phase (r ^ 0) from 
the decoupled phase (r — 0). Right: Schematic plot of the magnetization m-(h/Tf;). For Tq ^> 
(red dash dotted line), Tq ks T% (blue solid line), and Tq •< (green dashed line). The dotted line 
indicates the initial slope in the linear response regime, where m z (h) = h/Tg. The black dot refers 
to the complete polarization of the local Kondo spins, m z = 1/2, which occurs at the critical field 
h®. — Tk- Figures from ref. (4l . 

Ii/Tk- The low field regime is given by the linear response, m z {h) = h%i oc {T = 0), 
which, by definition, yields m z {li) — Ii/Tq. Since the critical field characterizing a 
full polarization of the local moments is of the order of the Kondo temperature, 
h® m Tk, we distinguish three typical cases. Tq w Tk is a standard situation, where 
m z increases linearly with h, until saturation. For To ^> Tk, the linear regime allows 
only a small magnetization when h < fP c . Thus, at about h as h° c , such systems exhibit 
a meta-magnetic transition from an unpolarized Fermi-liquid to the polarized spin 
lattice. In the opposite case, Tq <C Tk, the linear regime saturates around h ~ Tq 
h®. The intermediate regime Tq <h< Tk is expected to be non-universal. Eventually, 
a magnetization plateau can occur, at m z w (1 — n c )/2, similar to the one obtained 
for the entropy [see Fig . [2j . 

It is well known that the transition r = 0, defining Tk within the mean-field ap- 
proximation, becomes a crossover when more accurate methods are used. There are 
several examples of systems where a finite temperature crossover ends up to a quan- 
tum critical point, i.e. a transition at zero temperature. Wether the transition at the 
critical field h° c would survive beyond the mean-field is an open issue. 



4.2 Extra-RKKY interaction 

The possibility of destabilizing a Kondo phase in favor of a magnetically ordered 
ground state was first discussed by Doniach (8), by comparison of the Kondo tem- 
perature, T K ~ De-^^PoM, with the RKKY energy, J RKKY ~ p (/io)J|. 

The mean-field approximation for the Kondo lattice Hamiltonian ([T]i can hardly 
provide a correct description of a transition from a Kondo ground state to a phase 
with magnetic ordering. This results from the two following reasons: 
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( i) Difficulty in driving the system to a second order magnetic transition: this would 
require a mechanism leading to a continuous vanishing of the zero temperature ef- 
fective hybridization, r(T = 0). However, from the Kondo lattice Hamiltonian ([]]), 
the only possibility of vanishing r(T = 0) involves the non-interacting DOS, po(Mo) 
vanishing. This does not correspond to the physical situation of a magnetic transi- 
tion induced by the RKKY interaction. 

( ii) Difficulty in describing the criticality of transport properties: Within the mean- 
field approach, the system is either in a strong coupling regime, with r ^ and a 
heavy Fermi liquid ground state, or in a fully decoupled regime, with r = 0. Since 
the transport properties are governed by the conduction electrons, no universal non- 
Fermi-liquid behavior can be predicted at the mean-field level. 

These two difficulties can be understood from the large— formulation of the 
mean-field approximation: Coleman 11331 . and Read, Newns and Doniach [34] have 
shown that magnetic ordering involves at least processes of order 1 /N, i.e., fluctu- 
ations around the mean-field. Describing the criticality obtained from calculations 
taking into account the fluctuations and the emerging compact gauge field theory 
is beyond the scope of this presentation. We simply refer to the works of Senthil, 
Sachdev, and Vojta [40, 41], who analysed the possibility of a quantum phase tran- 
sition between a Fermi liquid phase with a large Fermi surface, and a (partially) 
fractionalized Fermi liquid phase. A Kondo breakdown quantum critical point has 
also been identified by Pepin B2l l43l . who obtained, for example, a specific ex- 
ponent diverging logarithmically in temperature, as observed in a number of heavy 
fermion metals. 

Whilst Difficulty (ii) can not be easily cured without the fluctuations, we de- 
scribe here how Difficulty ( i) can be fixed at the mean-field level. The Kondo lattice 
Hamiltonian ([U is generalized to a so-called Kondo-Heisenberg lattice, by adding a 
supplementary superexchange term, as follows: H H + L^/z/S/Sy. This general 
model has been introduced first by Sengupta and Georges ||44| . and studied later 
within various methods. In terms of auxiliary fermions, introduced in Section [2] 

the superexchange term is rewritten as JijSiSj = fiafia'f- a 'fjO' which de- 

scribes the spin-flip processes between two Kondo impurities, on sites i and j. Up 
to now, this description is very general and holds for any magnetic coupling Jfj, 
which can be periodic (ferro or antiferromagnetic), or randomly distributed (disor- 
der case). We will now describe two complementary mean-field approaches for the 
superexchange. 

The first one was introduced by Coqblin et al. B31 in the case of an antiferromag- 
netic nearest neighbor exchange 7y = Jaf < 0. In the paramagnetic Kondo phase, the 
intersite exchange is approximated within a Resonant Valence Bound decoupling: 
JijSiSj h-> rjj'Eaflcfjo, with Tij = Jj f- Za(fiaf]a)- The magnetic, Heisenberg-like, 
interaction generates an effective self-consistent dispersion for the /—fermions. Us- 
ing this method, Coqblin et al. show that the Kondo effect disappears abruptly for 
low band filling and/or strong intersite coupling [22. 23). 

We have developed another mean-field method 11461 . based on the DMFT, which 
we first applied to a disordered system where Jfj are given by a Gaussian distribution 
with an average [Jjj] = 0, and a variance [Jfj] =jj.ln this case, the superexchange 
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Fig. 9 (Color online) Schematic phase diagram of the disordered Kondo-Heisenberg model in the 
J c \ — T plane. Kondo temperature (red solid line) and coherence temperature (green dashed line) 
as function of Jj are shown for fixed values of J% and n c in the case of Tq < T^. The system is a 
heavy Fermi liquid below Tq(Jj). Above the line Tk(Jj), the localized spins are essentially free for 
Jd < T, whilst forming a highly correlated spin liquid for Jj>T '. All the lines represent crossovers. 
Figure from the analysis of Ref. |46| 

term generates a local energy-dependent self-energy for the /— fermions, which is 
determined by self-consistent relations similar to the ones obtained by Sachdev and 
Ye for a pure spin disordered system 1471 . We have also applied this method to a 
model with constant nearest-neighbor antiferromagnetic exchange (48), where the 
self-consistent equation for the /—fermions self-energy depends on the lattice struc- 
ture. In both cases, disordered or periodic model, we obtained a quantum critical 
point corresponding to the breakdown of Kondo effect, characterized by a vanish- 
ing of the effective hybridization, r(T = 0) = 0. Figure [9] depicts the phase diagram 
that we obtained for the disordered Kondo-Heisenberg model l46l . We also have ob- 
tained a similar phase diagram for the periodic, i.e., non-disordered, case f48l . Here, 
one relevant point is that the quantum critical point emerging from our mean-field 
approach does not necessarily correspond to the onset of magnetic ordering. This 
suggests an interpretation in terms of topological transition, without breaking of 
symmetry, but with a violation of Luttinger theorem, as discussed in Refs. I40ll4ll . 
Furthermore, 1 /N corrections can reveal a magnetic instability at a value of the cou- 
pling, J ( [ which is smaller than J c d = Tk . In this latter case, a symmetry breaking is 
expected. Whereas criticality in heavy fermions is governed by a Kondo breakdown 
critical point or by a magnetic transition, what is the nature of the non-Kondo phase, 
and whereas the Kondo breakdown coincides or not with magnetic ordering are still 
open questions. The answers probably depend on the system. 



Low energy scales of Kondo lattices: mean-field perspective 



15 



The quantum critical transition obtained here results from the competition be- 
tween the Kondo effect and the fluctuations of the /— fermions. These fluctuations 
are precisely the microscopical mechanism which is required to fix Difficulty (i), 
i.e., the difficulty in driving the system to a second order magnetic transition within 
the mean-field approximation. In the approach introduced by Coqblin et al, the 
/— fermion fluctuations are generated by an effective dispersion, which is connected 
to the Resonant Valence Bound decoupling. In our approach, the fluctuations are in- 
cluded within a self-consistent local self-energy. In both approaches, fluctuations 
are characterized by an energy, Jrkky = Jaf or Jj, and the Kondo effect disappears 
above the critical value Jrkky ~ Tk- 



5 Conclusions 

Important properties of Kondo systems can be obtained from the mean-field approx- 
imation. Some of them have not been presented here, like, for example, the effect 
of a pseudo-gap fl49ll50l . Also, for the sake of clarity, the description was restricted 
to the 'standard' Kondo model. Of course, the mean-field approximation has been 
generalized and applied to more realistic models, with, for example, a momentum- 
dependent hybridization between conduction electrons and / ions ||5T). 

Here we have focused on the low energy scales: the Kondo temperature, Tk, char- 
acterizing the temperature crossover below which conduction electrons and local 
moments are strongly coupled; the coherence energy, 7b, characterizing the Fermi 
liquid ground state; and Jrkky, the intersite magnetic correlation energy. At the 
mean-field level, Jrkky is negligible or neglected, and any study of magnetic crit- 
icality requires the inclusion of fluctuations around the mean-field. Nevertheless, 
the Kondo effect can break down at the mean-field level if the RKKY interaction is 
added 'by hand', generalizing the Kondo lattice to a Kondo-Heisenberg model. 

For the 'pure' Kondo lattice, the mean-field approximation provides explicit ex- 
pressions for Tq and Tk- Both quantities depend on the Kondo coupling with the 
same exponential factor, e _I /-^Po(Mo) f/his explains the strong sensitivity of these 
energy scales with respect to changes in the system (e.g. doping, or pressure). The 
ratio Tq/Tk, which does not depend on Jk, appears to be a promising quantity for 
analyzing universal behaviors of heavy fermion compounds. Whilst it is equal to 
one in dilute systems, Tq/Tk depends in fact on the electronic filling, n c , the band 
shape, and the impurity concentration, x. A universal regime with 7b <C Tk is ex- 
pected in the exhaustion limit, n c <C x, which might be difficult to access experi- 
mentally. More typical situations correspond either to the universal dilute regime, 
n c 3> x, or to the non-universal dense regime x ~ n c . In the latter case, the shape of 
the non-interacting DOS becomes relevant, and can lead to Tq <C Tk if the chemical 
potential, /1q, is close to a local maximum of the DOS, or to Tq ^> Tk if Ho is close 
to a local minimum. 

The experimental determination of Tq and Tk is straightforward for systems 
which are 'deeply' in a Fermi liquid phase: for example, one can determine Tq from 
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the specific heat Sommerfeld coefficient, and from the temperature dependence 
of the magnetic part of the entropy. The determination of Tg might become more 
tricky in the vicinity of a magnetic transition, where the freezing of the entropy is no 
longer due to Kondo singlet formation, but to magnetic intersite correlations instead. 
In this case, one should find another physical observable which would enable an un- 
ambiguous determination of Tk . In this case, a systematic experimental analysis of 
the ratio Tq/Tk might reveal interesting universal behavior. 

Some of the results obtained at the mean-field level have been confirmed by exact 
numerical methods. For example, DMFT calculations have shown that Tq/Tk does 
not depend on the Kondo coupling. One limitation of the mean-field is also well 
known: its weakness in describing criticality, where fluctuations become important. 
Nevertheless, some predictions of the mean-field have not been checked yet with 
more accurate methods (DMFT or 1 /N corrections). This is, for example, the case 
of the band shape effect. 

Acknowledgements I thank the organizers of the NATO Advanced Research Workshop on Prop- 
erties and Applications of Thermoelectric Materials. I acknowledge my collaborators on the works 
presented here: P. Fulde, A. Georges, M. Grilli, and V. Zlatic. I am also grateful to A. Rosch, M. 
Vojta, A. Klopper, H. Weber and N.B. Perkins for useful discussions regarding this manuscript, 
and to P. Nozieres, B. Coqblin, and C. Lacroix for fruitful advises. 



References 

1. P. Fulde, P. Thalmeier, and G. Zwicknagl, Strongly Correlated Electrons, Solid State Physics 
Vol. 60 (Elsevier, New York, 2006). 

2. H. v. Lbhneysen, A. Rosch, M. Vojta, and P. Wolfle, Rev. Mod. Phys. 79, 1015 (2007). 

3. A.C. Hewson, The Kondo Problem to Heavy Fermions (Cambridge University Press, Cam- 
bridge, England, 1993). 

4. S. Burdin, and V. Zlatic. larXiv;0812.1137l to appear in Phys. Rev. B. 

5. V. Zlatic, R. Monnier, J. Freericks, and K.W. Becker, Phys. Rev. B 76, 085122 (2007). 

6. V. Zlatic, R. Monnier, and J. Freericks, Phys. Rev. B 78, 0451 13 (2008). 

7. Y.F Yang, Z. Fisk, H.O. Lee, J.D. Thompson, and D. Pines, Nature 454, 611 (2008). 

8. S. Doniach, Physica B& C 91, 231 (1977). 

9. J.J. Joyce, A.J. Arko, J. Lawrence, PC. Canfield, Z. Fisk, R.J. Bartlett, and J.D. Thompson, 
Phys Rev. Lett. 68, 236 (1992). 

10. J.J. Joyce, A.J. arko, A.B. Andrews, and R.I.R. Blyth, Phys. Rev. Lett. 72, 1774 (1994). 

11. L.H. Tjeng, S.-J. Oh, C.T Chen, J.W. Allen, and D.L. Cox, Phys. Rev. Lett. 72, 1775 (1994). 

12. D. Malterre, M. Grioni, and Y. Baer, Adv. Phys. 45, 299 (1996). 

13. R. Mallik, E.V. Sampafhkumaran, PL. Paulose, J. Dumschat, and G. Wortmann, Phys. Rev. 
B 55, 3627 (1997). 

14. E.D. Mun, YS. Kwon, and M.H. Jung, Phys. Rev. B 67, 033103 (2003). 

15. J. M. Lawrence, P. S. Riseborough, C. H. Booth, J. L. Sarrao, J. D. Thompson, and R. Osborn 
Phys. Rev. B 63, 054427 (2001). 

16. E.D. Bauer, C.H. Booth, J.M. Lawrence, M.F. Hundley, J.L. Sarrao, J.D. Thompson, PS. 
Riseborough, and T. Ebihara, Phys. Rev. B 69, 125102 (2004). 

17. P, Nozieres, Ann. Phys. (Paris) 10, 19 (1985). 

18. P. Nozieres, Eur. Phys. J. B 6, 447 (1998). 

19. C. Lacroix, Solid State Commun. 54, 991 (1985). 



Low energy scales of Kondo lattices: mean-field perspective 



17 



20. C. Lacroix, J. Magn. Magn. Mat. 60, 145 (1986). 

21. S. Burdin, A. Georges, and D.R. Grempel, Phys. Rev. Lett. 85, 1048 (2000). 

22. B. Coqblin, M.A. Gusmao, J.R. Iglesias, C. Lacroix, A. Ruppenthal, and A.S.D. Simoes, 
Physica B 281, 50 (2000). 

23. B. Coqblin, C. Lacroix, M.A. Gusmao, and J.R. Iglesias, Phys. Rev. B 67, 064417 (2003). 

24. M. Jarrell, H. Akhlaghpour, and Th. Pruschke, Phys. Rev. Lett. 70, 1670 (1993). 

25. M. Jarrell, Phys. Rev. B 51, 7429 (1995). 

26. A. N. Tahvildar-Zadeh, M. Jarrell, and J. K. Freericks, Phys. Rev. B 55, R3332 (1997). 

27. A. N. Tahvildar-Zadeh, M. Jarrell, and J. K. Freericks, Phys. Rev. Lett. 80, 5168 (1998). 

28. A. N. Tahvildar-Zadeh, M. Jarrell, Th. Pruschke, and J. K. Freericks, Phys. Rev. B 60, 10782 
(1999). 

29. A. Georges, G. Kotliar, W. Krauth, and M.J. Rozenberg, Rev. Mod. Phys. 68, 13 (1996). 

30. W. Metzner, and D. Vollhardt, Phys. Rev. Lett. 62, 324 (1989). 

31. P. Nozieres, J. Phys. Soc. Jpn. 74, 4 (2005). 

32. C. Lacroix, and M. Cyrot, Phys. Rev. B 20, 1969 (1979). 

33. P. Coleman, Phys. Rev. B 28, 5255 (1983). 

34. N. Read, D.M. Newns, and S. Doniach, Phys. Rev. B 30, 3841 (1984). 

35. B. Coqblin, and J.R. Schrieffer, Phys. Rev. 185, 847 (1969). 

36. T.A. Costi, and N. Manini, J. Low Temp. Phys. 126, 835 (2002). 

37. S. Burdin, and P. Fulde, Phys. Rev. B 76, 104425 (2007). 

38. R.K. Kaul, and M. Vojta, Phys. Rev. B 75, 132407 (2007). 

39. J.A. Blackman, D.M. Esterling, and N.F. Berk, Phys. Rev. B 4, 2412 (1971). 

40. T. Senthil, S. Sachdev, and M. Vojta, Phys. Rev. Lett. 90, 216403 (2003). 

41. T. Senthil, M. Vojta, and S. Sachdev, Phys. Rev. B 69, 035111 (2004). 

42. C. Pepin, Phys. Rev. Lett. 98, 206401 (2007). 

43. I. Paul, C. Pepin, and M.R. Norman, Phys. Rev. Lett. 98, 026402 (2007). 

44. A.M. Sengupta, and A. Georges, Phys. Rev. B 52, 10295 (1995). 

45. J.R. Iglesias, C. Lacroix, and B. Coqblin, Phys. Rev. B 56, 1 1820 (1997). 

46. S. Burdin, D.R. Grempel, and A. Georges, Phys. Rev. B 66, 045 1 1 1 (2002). 

47. S. Sachdev, and J.W. Ye, Phys. Rev. Lett. 70, 3339 (1993). 

48. S. Burdin, M. Grilli, and D.R. Grempel, Phys. Rev. B 67, 121104 (2003). 

49. D. Withoff, and E. Fradkin, Phys. Rev. Lett. 64, 1835 (1990). 

50. L. Fritz, and M. Vojta, Phys. Rev. B 70, 214427 (2004). 

51. H. Weber, and M. Vojta, Phys. Rev. B 77, 125118 (2008). 



